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ABSTRACT 
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Certain geomagnetic storms exhibit, in addition to the usual 
initial sudden positive impulse, a subsequent sudden negative 
impulse . The former is normally ascribed to a shock wave in the 
interplanetary medium, and it has recently been suggested that the 
latter may be ascribed to a reverse shock convected away from the 
sun by the solar wind. 

If the velocity of efflux of gas from a source is supersonic 
(with respect to the source), if the velocity is instantaneously 
increased, and if certain subsidiary conditions are met, a pair of 
shock waves will be produced which propagate away from the source. 
The "fast" shock propagates away from the contact surface in the 
ambient gas (which was ejected from the source before the change 
in efflux velocity), while the "slow" shock propagates away from 
the contact surface in the driver gas, but has an outward velocity 
when this velocity is measured relative to the source . 

This problem (which may be identified with a classical problem 
considered by Riemann) is discussed in its relation to the produc- 
tion of pairs of shock waves by the enhanced solar wind produced 
by a solar flare . The equations giving the relationship between 
the velocities of the shock waves and of the ambient and driver 
gases become very simple in the strong-shock approximation. 
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It is shown that the propagation times of the positive and 
negative impulses of the July 11, 1959 magnetic storm may he 
explained satisfactorily on the basis of this theory. 

1 • INTRODUCTION 

Following the original suggestion of Gold [1955], it is now 
generally accepted that the positive sudden impulse which normally 
marks the beginning of a geomagnetic storm is caused by a shock wave 
in the interplanetary medium (now referred to as the "solar wind" 
[Parker, 1963] ) which, upon arrival at the earth, compresses the 
magnetosphere. Sonett and Colburn [1965] have recently drawn 
attention to the less frequent phenomenon of negative sudden impulse 
and suggested that this might be attributed to a "reverse" shock wave 
traveling away from the sun, relative to the sun itself, but toward 
the sun, relative to the solar wind. 

One of the possible mechanisms for the generation of such a 
pair of shock waves, in addition to other mechanisms described by 
Sonett and Colburn [1965], is that both shock waves originate near 
the sun as a result of a sudden increase in the velocity (and 
perhaps density) of the solar wind. This increase in flux is of 
course to be ascribed to the solar flare responsible for the geo- 
magnetic storm. Since the quiescent solar wind is supersonic 
(with respect to the sun), one expects that the enhanced solar wind 
produced by a flare will also be supersonic . The situation with 
which we are concerned is therefore, in its simplest representation. 


that of a driver gas being ejected from a source with supersonic 
speed (with respect to the source) into an ambient gas which has 
been ejected with supersonic speed from the same source . This 
situation already possesses one of the essential requirements for 
the production of a shock pair (consisting of a normal shock in the 
ambient gas and a propagating reverse shock in the driver gas ) , namely 
that the driver gas should have supersonic velocity with respect to 
the source from which the gas originates. The purpose of this article 
is to set up the simple basic equations for the phenomenon in order to 
determine what further conditions should be met, and to inquire into 
the likelihood of their being met 

This problem is closely related to a problem considered by 
Riemann and quoted by Courant and Friedrichs [ 19^8 (p. l8l)]: A 

straight tube of uniform cross section is divided into two parts by 
a plane diaphragm; to the left of the diaphragm is the gas we refer 
to as the "driver” gas and to the right of the diaphragm is the 
gas we refer to as the "ambient" gas . We now suppose that means 
exists for setting each gas instantaneously into uniform motion. 

Then the initial conditions on the problem are that the driver gas 
and ambient gas have pressure, density, and velocity p^, p^, v^ 
and p a , p a , v , respectively, and that, at the time these conditions 
are established, the diaphragm is removed. (See figure 1.) 

Riemann showed that the system will behave in one of four 
possible ways involving either a shock wave or a rarefaction wave 
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in each gas. However, for a rarefaction wave to be set up in each 
gas, it is clearly necessary that the gases recede from each other, 
that is, that v d < v a . Hence, if we stipulate that v d > v a , there 
must be either two shock waves (one in each gas) or a shock wave 
in one gas and a rarefaction wave in the other gas . 

Since we are concerned only with the possibility of two 
shock waves being established, this is the only case we shall discuss 
in detail. By considering that the strength of one shock wave is 
zero, we can find the condition which must be met in order that two 
shock waves should be set up. 

Figure 2 shows the situation in which a shock wave is excited 
in each gas. For later purposes, the shock wave in the ambient gas 
will be termed the "fast shock" and that in the driver gas the 
"slow shock." If we consider that v d > v a > 0, it is clear that 
the velocity of the contact surface, v c , must be positive. Hence 
Vf. s , the velocity of the fast shock wave, must be positive. However, 
the velocity v ss of the slow shock may have either sign. 

Hie case of particular interest to us is that v ss is positive. 
In this case, we may imagine that both the ambient gas and the 
driver gas are ejected by a source, with respect to 'which velocities 
are measured. Up to time t = 0, the pressure, density, and veloc- 
ity of the gas ejected by the source are p a , p a , v a ; for times 
subsequent to t = 0, the pressure, density, and velocity of the 
ejected gas are p , p d , v d - 


r 
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One would expect (and this is verified in the next section) 
that v ss is necessarily negative if the driver gas is in subsonic 
flow (and the conditions of the two gases are such that two shock 
waves are generated) . Hence a subsonic driver gas can never produce 
a shock pair, of the form shown in Figure 2, propagating away from 
the source . 

We now proceed to a discussion of the mathematical relations 
governing the production of a pair of shock waves by a supersonic 
driver gas . 


2 . THE SHOCK RELATIONS 


We now consider in more detail the configuration depicted in 
Figure 2 in which p^ is the pressure of the shocked driver gas, 
etc. The velocities of the shocked gases and of the contact surface 
are all equal. 


vl = v'-, = v 


( 2 . 1 ) 


a d ' c 

and the pressure is the same on both sides of the contact surface: 


Pa = ?d 


( 2 . 2 ) 


It is convenient to introduce Mach numbers in discussing the 


shock relations. We therefore introduce the speed of sound, c, where 

c = (7p/p) 1/2 (2.3) 


7 being the ratio of specific heats. We now characterize the 
velocities of the ambient gas and shocked ambient gas, with respect 
to the fast shock wave, by Mach numbers M a and M^_, writing 
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Similarly , we write 


fs " v a 

= M a c a 

(2.4) 

fs " v a 

= M a c I 

(2.5) 

d " v ss 

= M d c d 

(2.6) 

d - v ss 

= M!,c' 
d d 

(2.7) 


Since gas flows into a shock wave at supersonic velocity (measured 
relative to the shock wave) and emerges at subsonic velocity, M a 
and are greater than unity, whereas and are less than 
unity . 

We may note immediately from equation (2.6) that v < 0 if 

s s 

v d < c d> "that is, that the slow shock has negative velocity if the 
driver gas is in subsonic flow, a result referred to in 
Section 1. 

The following relations (the Rankine -Hugoniot equations 
[Courant and Friedrichs, 1948 (p. 129)]) hold at the fast shock. 


P a < v fs ‘ v a) = PA( v fs ' V ;L) 


( 2 . 8 ) 


P a ( v f s “ O* + Pa = P«( v - 


a' f s 


v^)“ + P ; 


f r* /■> \ 
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fs 


vJ 2 + 


■a 1 


7 - 1 P £ 


— ( v 
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fs 


vl ) 2 + — 1 — ^ 
a 7 - 1 P 


rr (2-10) 


and similar relations at the slow shock. 


The values of p, p, and M (or, alternatively, the values of 
p, c, and M) on one side of a shock wave determine the values on 
the other side. For instance, one may determine M* from M (or 
vice versa) from the relation 


m ,2 = (z - HI + 2 

2 7 M 2 - (7 - 1 ) 


( 2 . 11 ) 


Then the densities, pressures, 

£1 = (y + i)m 2 

P (7 - 1)M 2 + 2 


and sound speeds are related by 


( 2 . 12 ) 


p' _ 27M 2 - (7 - l) 
P 7+1 


( 2 . 13 ) 


£1 = {[27M 2 - (y - 1) ] [ (7 - 1)^ + 2 ] } 1/2 

c (7 + 1 )M 

If p, c, and v are given for the driver gas and for the ambient 
gas, the above relations for the contact surface and for the two 
shock waves provide nine equations for determining the nine 
unknowns p^, p^, c^, c£, v^, v*., v fg , v ss , and v Q . 

We may now inquire about the response of the magnetosphere to 
the arrival of a fact or slow shock wave . Quite detailed calcula- 
tions would be necessary for determining the precise change in the 
magnetosphere on the arrival of a shock wave. At this time, however, 
we wish simply to verify that there will be a compression of the 
magnetosphere, with a corresponding increase in the horizontal 
component of the magnetic field at the surface of the earth, at 
the time of arrival of the fast shock, and opposite changes at 
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the time of arrival of the slow shock. The calculations of Spreiter 
and Hyett [1963] show that the "side” and "tail” of the magnetosphere 
are sensitive to the value of p and insensitive to the value of 
pv 2 . Since p 1 > p, the side and rear of the magnetosphere should 
he compressed when the fast shock arrives and expanded when the 
slow shock arrives. 

One may expect further that the response of the front of the 
magnetosphere will he related to the total flux of momentum, which 
is given hy 

P = p + pv 2 (2.15) 


One nay verify that 

P* - P = (p« - p)v s 2 (2.16) 

■where v g is the velocity of the shock wave (v^ s or v ss ) . Since 
p' > p, P 1 > P. However, we may instead consider the stagnation 
pressure (the pressure at the "nose” of the magnetosphere) as an 
indicator of the influence of the solar wind on the front of the 
magnetosphere . For flow which is substantially supersonic with 
respect to the earth, we may use the formula [Landau and Lifschitz, 




I 


j 


p s = 



7- 1 2 

p-v 


(2.17) 


In this case, pv 2 is much larger than p, so that the change in 
pv 2 at the shock wave is given approximately hy (2.16) . The 
stagnation pressure will therefore increase with the arrival of 
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the fast shock and decrease with the arrival of the slov shock. 
These changes in pressure will lead to a sudden compression and 
expansion, respectively, of the magnetosphere, leading to a sudden 
increase and a sudden decrease, respectively, in the horizontal 
component of the geomagnetic field at the earth 1 s surface . 

We now consider the conditions which must he satisfied in 
order to obtain two shock waves rather than one shock wave and one 



Equation (2.19) may he expressed as a relation between 

c n/(vn - v ) and p /p . If we derive this relation and then note 
cL / x a a a' 

that higher relative gas velocities (or lower temperatures) clearly 
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favor the process of the formation of two shocks, we find that one 
of the conditions for the formation of two shocks is [Landau and 
Lifschitz, 1959 (p • 3^3)1 


v d " v a 


< 


-7 
2 ' 


(7 ♦ 1) + (7 - Djj 


'I 1/2 
>- 




- 1 


P 

if ^ > 1 (2.20) 


We see immediately that, if p /p < 1, the condition should 

3." CL 

he replaced by 


2 7 




V , - V 

& a 


< 


(7 + 1) ~ + (7 - 1) 

1 3 . 


- 1 


P, 

if p->l 


( 2 . 21 ) 


The condition represented by (2.20) and (2.21) is shown 
schematically in Figure 3 • It appears from this diagram that any 
significant sudden increase in velocity of the solar wind, if it 
is not associated with a remarkable change in pressure , is likely 
to produce the double -shock configuration. 


3 . STRONG -SHOCK APPROXIMATION 


The relations derived in the preceding section simplify if the 
Mach numbers M a and are large compared to unity. It appears 
that this requirement is met by the shock waves produced by the 
enhanced outflow of solar wind which occurs at the time of a major 
solar flare, such as produces a significant geomagnetic storm. In 
the following formulas, we adopt the value 7=5/3 appropriate to 
a fully ionized gas such as the solar wind. 




If M » 1 , we obt ain from (2.1l) 


/y n \ l/2 

M ' < 2 y ) =0 ' 45 

(3-1) 

Hence (2.12) leads to the familiar result 


P 7 - 1 

(3.2) 

( 2 . 13 ) leads to 

— « 27 M 2 - 1.25 M 2 

P 7 + 1 

(3-3) 

and ( 2 .lb) leads to 


M'c> 7-1 

— — a* = 0 .25 

Me 7+1 

(3-4) 

We now find from (2.4) and ( 2 . 5 ) that v., and v' 

is & 

may be 

related to v a , c , and M a by 

v f s = V a + M a c a 

(3-5) 

v a = v a + 0 -75 M a c a 

(3-6) 

Similar relations follow from ( 2 . 6 ) and (2.7). If, in 

addition, we 

use ( 2 . 1 ), these relations lead to 


v ss = v a + 0 - 75 M a c a " 0 - 2 5 M d c d 

(3-7) 

V d = v a + ° ,75 M a c a + 0 -75 M d c d 

(3-8) 


In order to find M a and M d , we need another relation in 
addition to (3-8). This is provided by the pressure condition (2.2). 
We write 



(3-9) 


is the ratio of the density of the driver gas to that 


so than A 1 2 

of the ambient gas. Equations (3-2) and (3.3) show that 

M a 2 P a = % 2 P d (3-10) 

from which, with the help of (2.3) and (3-9)^ we find that 


^a c a ~ ^% c d 


(3-11) 


We now find from (3*8) and (3*11) that 

EA 


M a c a 3(1 + A) 


(v d - 


(3-12) 


It is now possible to express the velocities of the fast and slow 

shock waves in terms of v^ and v^ as follows: 

a ci 


3 - A EA 

v h v - 


fs 3(1+ A) a 3(i + A) 


= 


v + 


3 ^ - 1 


v. 


uS 3(1 + ^) a 3(1 + A) d 


1 . DISCUSSION 


(3.13) 
(3 -14) 


It is now necessary to inquire whether the above formulas may 
be fitted to observations of geomagnetic storms displaying both 
positive and negative sudden impulses. It is necessary to find a 
geomagnetic storm for which the associated flare is well identified. 
For this reason, we select the flare of July 10, 1959 and the 
resulting geomagnetic storm of July 11, 1959- The positive sudden 


inpulse occurred 38.3 hours after the flash phase of the flare, and 
the negative sudden impulse occurred 8 .5 hours thereafter . 

When we come to relate our simple theory to observations of 
geomagnetic storms, some of the limitations of our treatment become 
evident. According to theory [Parker, 19^3], the speed of the 
solar wind does not vary greatly beyond a few solar radii from the 
sun. Although one would not expect that the enhanced flux of solar 
wind produced by an explosive phenomenon, such as a flare, would 
be constant over a long interval of time, the time scale of the main 
phase of a large solar flare has a duration of several hours, and 
this can be sufficient for discussion of the double -shock phenomenon. 
Since the Mach numbers of the two shock waves depend only on the 
relative velocity of the driver and ambient gases, on the density 
ratio of these gases, and on their temperatures, it is strictly 
necessary that all these quantities remain constant if the shock 
Mach numbers are to remain constant . One would expect the ratio 
of densities to remain constant, since each gas stream is expanding 
in the same way as the gas moves away from the sun. However, the 
temperatures of both gases will drop as the gases travel from the 
sun to the earth. For these reasons, and because we are neglecting 
the role of the magnetic field, one should not look for too close 
a correspondence between observations and the simple theories of 
Sections 2 and 3* It appears that correspondence will be better 
for strong shocks than for weak shocks, since the actual values of 
the temperatures of the ambient and driver gases do not appear 
in formulas (3*13) and ( 3 • 1^ ) • 


We therefore begin by applying equations (3-13) and (3*1*0 to 
the geomagnetic storm phenomenon. The above data concerning the 
time delay between the occurrence of a flare and the occurrence of 
the sudden impulses lead to uhe estimates Vf S = 103p km/sec and 
v ss = 888 km/sec. We readily find from (3*13) and ( 3 - 1*+ ) that 


v a = 3(v fs - v RR ) 


ss' 


(*v.l) 


so that, in our case, v. - = 591 km/sec, and that 


A = 


V„ - Y 

is a 


3v f s ■ W ss + v a 


(4.2) 


If v 0 is known, v, and A are now determined. Since no measure- 

a u. 

ments of the solar wind were made in 1959 > we must consider various 
possible values of v . In this way, we find that the observed 
time delays may be explained by any of the sets of parameters 
listed in Table 1. The Mach numbers given in Table 1 have been 
estimated on the assumption that c a = c a = km/ sec, corresponding 

to a temperature of 10 5 deg of fully ionized hydrogen gas . It is 
clear that, with this choice of temperature, the strong-shock 
assumption is unacceptable for v a = 400, 500 km/ sec , but is a fair 
approximation for v a = 600 km/ sec . A smaller value of the gas 
temperature would, of course, make the strong-shock assumption 
more acceptable . 

It is possible to draw up a similar table on the basis of 


the formulas of Section 2, provided one can relate Me to M'c' 
for each shock. We find from (2.11) and (2.1*0 that 



for 7 = 5/ 3- The procedure is as follows: 

For known values of v, , v and assumed values of v a , c 0 , c„ 

is ss a- 7 a 

one may determine Ma from (2.1) and then, using (4.3), one may 
determine v' (which is the same as v and v ' ) from (2.5). One 

^ ^ CL 

may then determine M^c^ from (2.7) and, determining from (4-3) 

one may finally find v^ from (2.6). If uhis procedure is applied 
to the event considered above in the strong-shock approximation, 
we arrive at the sets of parameters given in Table 2. The density 
ratio is evaluated from 

P d 27M a g - (7 - 1) e/ 5M a 2 - 1 c a 2 

P a 271^ 2 - (7 - 1) c a 2 5% 2 - 1 c d 2 

which is derivable from (2.2), ( 2 . 3 ), and (2.13). We note that 
for v a = 600 km/ sec, the estimate of the driver gas velocity is 
not very different from that obtained in Table 1, although there is 
a substantial difference in the estimate of P^/p a - For 
v o = 500 km/sec, there is an ambiguity in determining M^. For 

v a = 400 km/sec, no solution is possible. 

One may verify from Figure 3 , that, for all sets of parameters 
listed in Tables 1 and 2, the criterion for the generation of two 
shocks is satisfied. 

Although many approximations have been made in the course of 
this analysis, it appears from this discussion that positive and 


negative sudden impulses may, in fact, be attributable to a pair 
of shock waves generated by an enhanced flux of the solar wind 
at the time of a solar flare • 
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TABLE 2. Possible Parameters for July 11, 1959 Event 
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Fig. 1. 
Fig. 2. 
Fig. 3* 


Initial configuration in Riema nn’s problem. 
Configuration after collision of gas streams . 
Condition to be satisfied to oroduce two shock waves 
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